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Methods of extraction of the symmetry energy (or enthalpy) coefficient to temperature ratio from 
isobaric and isotopic yields of fragments produced in Fermi-energy heavy-ion collisions are discussed. 
We show that the methods are consistent when the hot fragmenting source is well characterized and 
its excitation energy and isotopic composition are properly taken into account. The results are 
independent of the mass number of the detected fragments, which suggests that their fate is decided 
very early in the reaction. 

PACS numbers: 21.65.Ef, 24.10.Pa, 05.70.Fh, 25.70.Mn 



I. INTRODUCTION 

In the past few years the importance of the symmetry 
energy term in the nuclear equation of state has stimu- 
lated a growing interest in isospin effects in nuclear reac- 
tions. Understanding the properties of asymmetric nu- 
clear matter both at normal densities and at densities 
away from the saturation density has an important im- 
pact on the study of the nuclear structure close to the 
drip lines 0] and on the study of astrophysical processes 
01- 

Recent measurements of the giant dipole [3fl, Pygmy 
dipole [4] and giant monopole [5] resonances in neutron- 
rich nuclei, neutron/proton emission 0], isospin diffusion 
01 and fragment isotopic ratio 0, 0] have provided initial 
constraints on the density dependence of the symmetry 
energy at sub-saturation densities. Refinement of these 
measurements with both stable and rare isotope beams in 
the near future will provide further stringent constraints. 
New crucial experimental constraints on the symmetry 
energy at supra-saturation densities are expected from 
the recent measurement on neutron-proton elliptic flow 
performed at GSI [10]. Nevertheless the conclusions are 
model dependent, different analysis methods provide dif- 
ferent results and a consistent description of all avail- 
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able experimental data employing one type of equation 
of state is still lacking. 

In this work, among the experimental observables com- 
monly used to explore the symmetry energy at sub- 
saturation densities, we will concentrate on yield ratios 
of fragment produced in multifragmentation processes at 
Fermi energies. Different methods to extract the "sym- 
metry energy" coefficient from this observable have been 
proposed. Among those, we will focus on three seem- 
ingly different approaches: isoscali ng fllT - il6| . m-scaling 
[171 ] and isobaric yield ratio method [18] . and we will show 
that consistent experimental results can be extracted pro- 
vided that the properties of the hot fragmenting source, 
such as its isospin composition and excitation energy, are 
properly taken into account. The paper is structured 
as follow: section IlTI briefly describes the three methods. 
Sections llIII and llVl present the motivation that led to this 
work and the experimental apparatus, respectively. The 
importance of the source reconstruction and the results 
are presented in sections [V] and IVII respectively, while 
section IVHI concludes. 



II. THEORETICAL BACKGROUND 

The study of the multifragmentation process in vio- 
lent heavy-ion collisions at Fermi energies is important 
for the investigation of the symmetry energy. During the 
multifragmentation process, which has been related to a 
nuclear phase transition |l9j-22]. subsaturation density 
may be achieved [23| . Free energies play an important 
role in mixed phase environments. Indeed the isotopic 
distribution of fragments produced in such collisions is 
governed by the free energy at the pressure and temper- 
ature of the fragmenting source. Thus, assuming that 
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the fragment production is governed by a purely statis- 
tical process at constant temperature T and pressure P, 
the yield of a fragment, with N neutrons and Z protons, 
Y(N, Z,T, P), can be related to the nuclear Gibbs free 
energy G(N,Z,T,P) pj|: 



Y a A~ 



exp 



Y(N,Z,T,P) = 

|_ G(iV,Z,T,F) + fe iv + ^ z | 



(1) 



where Yq is a constant, \i n and fj, p are the neutron and 
proton chemical potentials, T is the temperature of the 
emitting source. The factor A~ T originates from the en- 
tropy of the fragment, i.e. the Fisher entropy [24l . [25J. 
The use of the Gibbs free energy [26] depends on the 
validity of the assumption of fragment production via 
an equilibrium mechanism at constant pressure. Other 
assumptions, and therefore other thermodynamic state 
functions, could be used [23, [28j], but it is undoubtedly 
true that the true fragment production scheme requires 
a kinetic treatment. For instance if the volume is kept 
constant (freeze-out hypothesis) , then the Helmholtz free 
energy, F(N, Z,T,V), should be used. This ambiguity 
casts some doubts on the derived quantity, i.e. symme- 
try energy or enthalpy. 

Within a liquid drop description, the nuclear free en- 
ergy, F(N, Z, T, V), can be parametrized as a sum of the 
bulk, surface, Coulomb and symmetry free energy contri- 
butions. The symmetry energy term is usually expressed 
asflllll 



E sym {N, Z,T,V) 



C sym (T, V) 



(N-Z) 2 
A 



(2) 



where C syrn (T,V) is the symmetry energy coefficient, 
V is the volume and A = N + Z . Similarly, starting 
from the Gibbs free energy, G(N,Z,T,P), the sym- 
metry enthalpy term could be espressed introducing 
a symmetry enthalpy coefficient Chsym, depending on 
(T, P). The logic employed in this work, which assumes 
an equilibrium process at constant pressure fl4j . will 
lead to the extraction of the symmetry enthalpy. The 
assumption of a constant volume [l2| would have lead 
to the extraction of the symmetry energy. Experimen- 
tally whether the equilibrium process takes place at 
constant pressure or volume (freeze-out hypothesis) 
is not determined, therefore the ambiguity on the 
extracted quantity is retained. However some estimates 
show that the difference between the two quantities 
should be small below the critical point of the phase 
transition [30]. Above it, the PV product becomes 
finite, and the two quantities might significantly differ. 
Nevertheless in Ref. [27j it is shown that the canonical 
and grand-canonical ensembles predict similar results if 
one is only interested in the ratio of population of two 
adjacent isotopes, on which we will focus in this work. 
Keeping in mind this ambiguity, from now on we will 
refer to the experimentally extracted quantity as the 
"symmetry energy". 



We now give a brief review of the three methods that 
we will use in our experimental analysis. 
It has often been experimentally observed that the ratio 
of the yields of a fragment with N neutrons and Z pro- 
tons produced in two similar reaction systems with dif- 
ferent neutron-to-proton ratios is exponential in N and 
Z [Til - [l6j . In the grand canonical approximation, using 
Eq. Q] this ratio can be written as 



R 21 (N, Z,T, P) 



Y 2 (N,Z,T,P) 
Yi(N,Z,T,P) 



= Cexp{aN + (3Z) 

(3) 

assuming that the thermodynamic state points of the 
two equilibrated sources in the two reactions are the 
same. (The indices 1 and 2 denote the neutron-poor and 
neutron- rich system, respectively). This relation, known 
as isoscaling, has been found to describe the measured 
ratios over a wide range of complex fragments and 
light particles rather well [l2l | and to be a phenomenon 
common to many different types of heavy-ion reactions 
d, [H, GJ, Hi EH. This suggests that free energy 
components sensitive to the neutron-proton concentra- 
tion differences play a key role in the fragment formation 
process. 

Recently, the Modified Fisher Model of Ref. [H] has 
been used to interpret multifragmentation data and ex- 
tract information on the symmetry energy 
In this model, the free energy per particle is modified as 

nm f ,A,T,H)« G ( N > Z > T >V-^ N -^ Z (4) 

where rrif (— N ^ Z ) is the relative isospin asymmetry of 
the fragment. This is the Landau free energy. The ansatz 
here is that near a critical point all the dependencies of 
the free energy are contained in the order parameter m / 
and its conjugate field H [13, HH ■ In principle we do 
not need to specify if the system is at constant volume 
or pressure, as before. Comparing to Eq. [T] we should 
notice that ^>(rrif, A, T, H) includes the neutron and pro- 
ton chemical potentials. Within this approach, the ratio 
of the free energy per particle to the temperature near 
the critical point is described by the expansion [33T - [35j : 



^{mj,A,T,H) 



f + 0(m 8 f ). 



(5) 



The parameters a, b and c, which depend on T and p 
[IH, are used for fitting. We stress here that the "exter- 
nal field" in our case might be attributed to the difference 
in neutron and proton chemical potentials of the source 
[SHU- Indeed in the limit where fi p = —fi n , combining 
Eqs. S] and [51 we obtain 2H — /i„ — fi p . Other terms, 
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such as the volume, surface, Coulomb and pairing con- 
tributions, might be negligible near the phase transition, 
as suggested by experimental data (33| . 

An exponential dependence on has been experi- 
mentally observed for the ratio, i? 2 i, of the yields of 
a fragment (to/, A) produced in two similar reactions 
with different neutron-to-proton ratios (l7| . This scaling 
has been referred to as m-scaling, to distinguish it from 
the known isoscaling. Indeed the fragment yield ratio 
between two systems at the same thermodynamic state 
point depends only on the external field H/T: 



R 21 (m f ,A,T, H) 



A H 

Cexp( — m f A) (6) 



where AH/T = H 2 /T - Hi/T. Comparing Eq. Hand© 
and assuming that a = — 8 , we obtain = a, as has 
been suggested in Ref. [17]. As w will show that this 
relation is approximately satisfied in our experimental 
data, ignoring the isospin symmetry breaking Coulomb 
effects is a reasonable approximation for light nuclei [33| . 

Following the statistical interpretation of the isoscal- 
ing with SMM [l3[ and with the expanding emitting 
source model [l2| , it has been proposed that the symme- 
try energy coefficient can be extracted from the measured 
isoscaling parameters through the approximate formula 
El: 



C i sym (-^) 



4A' 



(7) 



The two models of Refs. [H [lj] are based on different 
assumptions, therefore we restrain from stating explicitly 
the C sym dependence on V or P, which is model depen- 
dent. 

In the literature a variety of definitions of the quantity 
A has been proposed [HI, [3. Among those we will focus 
on two commonly used expressions and on one recently 
suggested. 

In the context of SMM and the expanding emitting 
source model, A is the difference of the asymmetries of 
the equilibrated emitting sources, defined as [l2j: 



^source 



(8) 



In general, due to the phenomenon of fractionation, 
fragments may not have the same isospin ratio of the 
fragmenting source and this could affect the value of C sym 
obtained from the isoscaling coefficient. 

A study of the isospin fractionation and the isoscaling 
with AMD simulations [36] has shown a linear relation 
between a and (Z/A) 2 of fragments. The quantity A has 
then been expressed as 



A 
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where Z/(A) is the proton fraction of the most probable 
isotope for fragments of a given Z with A > 4 ("liquid 
part") [Til. 

More recently, R. Tripathi et al. [13] have shown, in 
an investigation of the nuclear phase transition using the 
Landau Free energy approach, that the position of the 
central minimum of the free energy is related to the av- 
erage isospin asymmetry of the fragments produced in 
each event, excluding neutrons and protons, rnj. This 
suggests A may be written as 
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(10) 



where (mf) is the event average of the fragment relative 
isospin asymmetry. 
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where N is the number of events and K the event multi- 
plicity excluding neutrons and protons. 

We would like to stress here that the definitions of 
A proposed in Eq. [8j [9] and [10] differ significantly. 
While Eq. [H] relates A to the characteristics of the two 
sources, Eqs. [9] and [10] express A as a function of the 
characteristics of the fragments. Moreover Eq. [101 as 
opposed to Eq. [5] takes into account, event-by-event, 
the multiplicity in addition to the composition of the 
fragments. Since the determination of C syrn /T from 
scaling parameters depend on A, the correct determina- 
tion of A is critical. 

Recently, a different method to determine C sym /T has 
been proposed, which does not depend on the definition 
of A. Within the Modified Fisher Model 0] it has been 
shown that the isotope yield ratio between two isobars 
differing by 2 units in the neutron excess I = N — Z and 
produced by the same source can be written as [l8j : 



~^S m i X (I + 2, A) — S m i X (I , A)} 



f rp ^ W(I+2,A,T,p)-W(J,A,T,p) + (p, n -n p ) 



(12) 



where W(I, A, T, p) is the free energy of the cluster at 
temperature T, which in Ref. [181 ] was approximated 
by the generalized Weiszacker-Bethe semi-classical mass 
formula [38l [39| . and S m i X {I,A) is the cluster mixing 
entropy. Apart the change of the notation for the free 
energy, we stress that a similar formula could be derived 
at constant pressure, thus the ambiguity discussed above 
remains also in this case. Note that the variable /, intro- 
duced for consistency with Ref. [IS], could be expressed 
as a function of the isospin asymmetry m f as / = m f A. 
The authors show that the symmetry energy coefficient 
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to temperature ratio can be expressed as: 

£^21^-^[lnR(3,l,A)- 

-lnR(l,-l,A)-5{3,l,A)} 

where 6(3,1, A) is the difference in the mixing entropies 
and can be neglected, being rather small compared to the 
other terms in Eq. [T3] [lj] . 

Equation [JjJ] can be derived within a Landau free en- 
ergy approach, which is used here to show that the depen- 
dence of C sym /T on the source characteristics, contained 
in H/T, cancels out to first order. The yield of each frag- 
ment can be expressed as a function of its free energy per 
particle F, as in Eq. [5] Neglecting 0(m 4 ) and choosing 
1 = 1, —1, i.e. mi = 1/A and m,2 = —1/A, and I = 3, 1, 
i.e. mi = 3/ A and mi = 1/A, the logarithm of the two 
yield ratio can be written as: 



lnR(l,~l,A) s= 
lnR(3, 1, A) w -A ri 



o A H 1 
ZA TA 



-a-§- 

2 U A 2 



(14) 



the experimental conditions, the effect of the secondary 
decay and the method of determination of A. 

In a previous work, S. Wuenschel et al.[40j have 
presented data from quasi-projectile fragmentation in 
78,86 Kr+ 58,64 Ni at 35A MeV, showing an independence of 
the isoscaling-extracted C sym /T values from the charge 
of the analysed fragments. Stringent constraints were ap- 
plied for a careful characterization of the source with a 
well defined neutron (proton) concentration. The data 
were collected by a 4tt charged- particle detector array. 
In the works presented in [17], lla, [29J, data for the re- 
actions 64 ' 70 Zn, 64 Ni+ 58 < 64 Ni, 112 ' 124 Sn, 197 Au, 232 Th at 
40 AMeV analysed with isoscaling, m-scaling and isobaric 
yield ratio methods instead show a strong dependence 
of C S y m /T values on the charge of the analysed frag- 
ment. The data were collected with a detector telescope 
placed at 20°. The telescope allowed the clear identi- 
fication of typically 6 — 8 isotopes for atomic number, 
Z, up to Z = 18. In contrast to the experimental con- 
ditions of S. Wuenschel et al. [4(|, the information on 
the fragmenting source was not available and the angu- 
lar coverage was limited (15° < 9 < 25°). Thus the yield 
of each isotope was evaluated using a moving source fit, 
selecting fragments produced in the fragmentation of the 
overlapping region of the two colliding nuclei. 

It should be mentioned here that, if compression effects 
and/or higher order terms in Eq. [TBJare important, the 
isoscaling and the isobaric yield ratio methods might pro- 
duce different results. Moreover, a different beam energy 
and different colliding nuclei could access different densi- 
ties thereby providing complementary information. Thus 
investigations under precise experimental constraints for 
the source properties are necessary, which is one of the 
goals of this paper. 

Furthermore, the observation of an increasing behav- 
ior of C S y m /T as a function of the fragment mass raises 
the question of the effects of the secondary decay process 
on the observables. Indeed, in previous works, excitation 
energies of the primary fragments estimated from the as- 
sociated light charged-particle multiplicities fill . IHJ sug- 
gest that effects of secondary evaporation could be im- 
portant. Another possibility could be that we are really 
measuring an enthalpy rather than the internal symmetry 
energy and different nuclei could be emitted at different 
volumes. 

The third issue is the expression of the quantity A 
that is used to extract C sy m/T from the isoscaling and 
m-scaling parameters. Indeed, while A from Eq. [5] has 
been widely used after its introduction [12], the study 
of the isospin fractionation suggests that the connection 
III. MOTIVATION OF THE PRESENT WORK between C sym /T and the scaling parameters should be 

found in the fragment isotopic asymmetry (Eqs. |TJ] and 
ITU]) . Whether the event composition (Eq. ITU]! or the 
average composition of each fragment (Eq. [TJ]) is more 
relevant is still an open question, which is specifically 
addressed in this work. 

In this work we investigate the importance of the 
source reconstruction and the calculation of A by com- 



where l/2a, fitting parameter defined in Eq. [5] is by 
definition C sym /T IM HE 13- Th e difference be- 
tween the two relations, as in Eq. [T_3l cancels out the 
dependence on ^ , removing the dependence on the char- 
acteristics of the source. 

We would like to stress here that the main difference 
between the first two methods of extracting C sym /T 
and the third is that in the first two cases we are com- 
puting the yield ratios of the same fragment produced 
by two different sources, while in the third case we 
consider different isobars produced by a common source. 
Therefore, in the first two cases the determination of 

C sym 

/T requires an estimation of while in the third 
case the dependence on the source characteristics is 
removed by computing the difference of two yield ratios. 
However, the isobaric yield ratio method is dependent on 
the selected isobars and vulnerable to secondary decay 
effects. Isoscaling and m-scaling are less affected by 
secondary de-excitation effects, since they are removed, 
to first order, when taking the ratio of the yields of the 
same fragment produced in two sources. We will show 
that the three methods, which are consistent within the 
Landau free energy approach (l8l . [29j . give C sy m/T in 
good agreement. Moreover the agreement of the results 
could suggest the method of determination of A. 



In the present work, fragment yield data from quasi- 
projectile fragmentation in 64 Zn+ 64 Zn, 70 Zn+ 70 Zn and 
64 Ni+ 64 Ni at 35A MeV have been analysed using the 
isoscaling, the m-scaling and the isobaric yield ratio 
methods. There are three issues for the determination 
of the Csym/T values in Eq. [71 and [T31 the importance of 
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Mass Identification Parameter |-p 



Figure 1. (Color online) Isotopic distributions for Z = 4 (left) 
and Z = 14 (right) for a Super- Telescope. The x axis is 
obtained by linearizing AE — E spectra as described in 
The gaussian fits (black, dashed) show the overlapping of the 
mass distributions for different isotopes. The overall gaussian 
fit is also showed (red, solid) for each Z. 

paring the C sym /T values obtained with the three differ- 
ent methods. 



IV. THE EXPERIMENT 

The experiment was performed at Texas A&M Univer- 
sity K500 superconducting Cyclotron. Beams of 64 Zn, 
70 Zn and 64 Ni at 35AMeV beam energy were impinged 
on 64 Zr^ 70 Zn and 64 Ni targets. The 4tt NIMROD-ISiS 
array [4J , |45| was used for the detection of charged par- 
ticles, obtaining isotopic resolution for Z < 17. The de- 
tector telescopes, arranged on 15 rings centered on the 
beam axis, were composed of one silicon detector backed 
by a CsI(Tl) crystal with PMT readout. Two telescopes 
per ring, referred to as Super- Telescopes, made use of 
two silicon wafers in front of the CsI(Tl), in order to im- 
prove the isotopic resolution. The charged-particle array 
was housed inside the TAMU Neutron Ball [45[, which 
measured the free neutron multiplicity. 

The mass identification of fragments is essential for 
isospin physics analysis. For Z < 2, mass identifica- 
tion was performed by pulse shape discrimination of the 
CsI(Tl) signal. The AE - E method was used for Z > 3 
on Si-Si and Si-Csl telescopes. In Fig. [T]the isotopic dis- 
tributions of Be and Si are shown as examples. We can 
see that the global fit (red, solid) reproduced the over- 
all spectrum well. The individual gaussian fits (black, 
dashed) show the overlapping of the mass distributions 
for different isotopes. Particles were confidently assigned 
a mass if the contamination from neighbour isotopes was 
less than 5 — 10%. As evidenced by the figure, the mass 
determination efficiency decreases for large A near the 
upper limit of the electronics range. Further details on 



Figure 2. (Color Online) Reconstructed quasi-projectile m s 
distribution from the 64 Zn+ 64 Zn (squares), 64 Ni+ 64 Ni (dia- 
monds), 70 Zn+ 70 Zn (triangles) reactions. The total yield of 
each system was normalized to unity. The m s distribution 
of the quasi-projectile obtained from the combined systems 
is plotted as stars. The total yield of this distribution was 
normalized to 3. 



the mass identification procedures and on the experiment 
can be found in (43l | . 

Particle and event selections were performed to se- 
lect quasi-projectile fragmentation events. The quasi- 
projectile source was reconstructed accepting in each 
event fragments with a longitudinal velocity relative to 
the largest fragment within the range ±65%, 60% and 
45% for Z = 1, Z = 2 and Z > 3, respectively [H]. This 
cut, later referred to as Vcut, is intended to remove frag- 
ments from non-pro jectile-like sources. Furthermore the 
total Z of the detected fragments included in the recon- 
struction was constrained to be in the range Z = 25 — 30 
(SumZ). Finally limits were placed on the deformation of 
the source, as measured by the quadrupole momentum, 
to select a class of events that are, on average, spher- 
ical. The quadrupole momentum, calculated from the 
measured particle momenta in the quasi-projectile frame, 
X^Pfj/EjPIjj was required to be less than 2 (Qcut). 
Details on the source reconstructions can be found in 
Ref. j4Z[. Free neutrons measured by the Neutron Ball 
were used to correct for the free neutrons emitted by 
the quasi-projectile using the procedure discussed in ref. 

V. IMPACT OF SOURCE RECONSTRUCTION 

Figure [2] shows the distribution of the reconstructed 
quasi-projectile isospin asymmetry m s — Na ~ Zs for each 
of the three systems and for the three systems combined 
together. In the m s expression, N s , Z s and A s are, re- 
spectively, the number of neutrons, protons and total nu- 
cleons in the reconstructed quasi-projectile. As expected, 
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Figure 3. (Color Online) Isoscaling parameter a extracted for 
Ami = 0.185 and 3.5 < E* t /A < 5 AMeV from each reaction 
and for the combined systems as a function of the fragment 
mass (A = 2Z). 



the systems 70 Zn+ 70 Zn and 64 Ni+ 64 Ni, which have very 
similar N/Z values (1.33 and 1.29 respectively), show al- 
most overlapping m s distributions, with average values 
mj = 0.15 and 0.13 respectively. The 64 Zn+ 64 Zn m s dis- 
tribution is shifted toward lower m s values (mj = 0.09), 
since the system is less neutron rich (N/Z = 1.13). The 
m s distribution of the three combined systems is cen- 
tred around m s = 0.12. The widths of the distributions 
are large compared to the difference in the average m s 
between the reacting systems. In Ref. [4Cj it has been 
shown that an improved isoscaling can be obtained when 
selecting narrow bins in N/Z of the fragmenting source, 
rather than performing a system-to-system isoscaling. In 
analogy, and extending that work, we performed both 
isoscaling and m-scaling between two different m s bins. 
For the isoscaling analysis using m s , AmJ refers to the 
difference of the mean m s values in the two m s bins. To 
allow a comparison of the results from the three different 
methods, for each m s bin combination the isobaric yield 
ratio was computed for fragments produced by a source 
with (m s ) = AmJ ± 0.0625. 

Fragment yield data was divided according to the rel- 
ative isospin asymmetry, m s , of the reconstructed source 
in 4 bins: m Sa = -0.0375±0.0625, m Sb = 0.0875±0.0625, 
m Sc = 0.2125 ±0.0625 and m Sd = 0.3375 ±0.0625, whose 
boundaries are also plotted in Fig. [2j The transverse 
excitation energy of the reconstructed source (El/) was 
calculated through calorimetry, as E3| : 



M C p 

E 



K? p (i)+M n (K?)-Q. 



(15) 



The first term is the sum of the transverse kinetic en- 
ergies in the quasi-projectile center-of-mass of the parti- 
cles (Mcp) belonging to the quasi-projectile (see SecJIVI 
- VCut). The free neutrons contribution to the excita- 
tion energy was estimated as the average neutron kinetic 
energy (if™) multiplied by the neutron multiplicity M n . 
The average transverse kinetic energy of the neutrons was 



calculated as the proton average transverse kinetic energy 
corrected for the Coulomb barrier energy [48]. The last 
term in the equation is the reaction Q value. The mass 
of the quasi-projectile was calculated as the sum of the 
masses of the charged particles belonging to the source 
and the neutron multiplicity. The transverse excitation 
energy per nucleon (E* /A) of the reconstructed source 
was calculated and the data was divided into five bins: 
(1.75 ± 1.75)MeV, (4.25 ± 0.75)MeV, (5.75 ± 0.75)MeV, 
(7.25 ± 0.75)MeV and E% /A > 8MeV. 

Figure [3] shows the isoscaling parameter a as a func- 
tion of the mass fragment A — 2Z extracted from each 
reaction (open symbols) and for the combined systems 
(full circles). The error bars are comparable to the 
size of the points. The isoscaling was performed be- 
tween two sources with different m s (m Sl = —0.002 and 
mJ7 = 0.183) so that AmJ = 0.185. The source excita- 
tion energy was required to be between 3.5 and 5AMeV. 
We observe that consistent a values for each mass were 
obtained for each system and for all the systems when the 
data was gated on m s . The same results were obtained 
for all the m s bin combinations and for all the excitation 
energies. Moreover a similar behavior was observed for 
AH/T extracted from the m-scaling. Thus the systems 
were combined together to increase the statistics. 

Significantly better isoscaling and m-scaling were 
obtained once two sources with different m s rather than 
two different reactions were selected. In Fig. 2] the 
natural logarithm of the yield ratio of Eq. [3] and -jlnR 2 i 
of Eq. UJ] are plotted as a function of N and m / , respec- 
tively. On the left the ratios are computed between the 
yields of fragments produced in two different reactions 
( 64 Zn+ 64 Zn and 70 Zn+ 70 Zn), while on the right the 
fragments are produced by two sources with different 
m s (mJT = 0.097 and rn^ = 0.180), to construct the 
neutron-poor and the neutron-rich sources, as required 
from Eqs. [3] and [5] Comparison of the left and the right 
panels demonstrates the improvement of isoscaling and 
m-scaling with a narrowly defined m s source. Indeed 
the isotopes line up on parallel, equally spaced lines in 
the isoscaling and on a single line in the m-scaling, as 
predicted by Eqs. [3] and [51 respectively. 

To investigate the impact of the source reconstruction 
on the extracted values, the source constraints described 
in this and the previous sections have been imposed se- 
quentially, on all the data, in order to illustrate the effect 
of each cut individually. The three methods have been 
applied to the so-selected data. 

Figure \5\ shows s J m , a and AH/T extracted from the 
isobaric yield ratio (Eq. [13)) . the isoscaling (Eq. [3]) and 
the m-scaling (Eq. [6]) methods, respectively, as a func- 
tion of the fragment mass number A. For the x-axis of 
the isoscaling and m-scaling methods A = 2Z. The a 
and AH/T parameters are determined both by individ- 
ual fits to the yield ratios for isotopes with a given Z. 
We first focus on the values obtained comparing 
70 Zn+ 70 Zn to 64 Zn+ 64 Zn systems (labelled as SysToSys 
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in figure) , to analyze the impact of a limited angular cov- 
erage and of the quasi-projectile reconstruction. We will 
refer to these values as system-to-system values. A gener- 
ally flat behavior is observed for a and AH/T, indepen- 
dent of the fragment mass number A, when no constraints 
are imposed on the data (stars). The trends are modi- 
fied neither by restricting our analysis to a limited an- 
gular range (diamonds) , nor by selecting quasi-projectile 
fragmentation events (i.e. applying the Vcut, SumZ and 
QCut - circles, triangles and full squares, respectively). 
This suggests that the quasi-projectile selection does not 
significantly bias our determination of a or AH/T. 

In contrast to the flat behavior of the isoscaling and 



m-scaling, 



increases from 9 to 20 over the A range 



of 11 to 31, when no constraints are imposed on the data 
(stars). The increase is slightly amplified for A = 25 and 
27 by selecting a limited angular range (diamonds) to 
simulate the angular coverage of Ref . [18j . The observed 
increase of 11 units in is similar to the one of 8 

units reported in Ref. [18| . Data for masses A > 27 in a 
limited angular range are not available due to a statistics 
limitation. We should notice that the angular selection 
does not affect the values extracted from the isoscaling 
and m-scaling, since we take the ratio of the yields of 
the same fragment (Eqs. [3] and [6]). This is not true for 
the isobaric yield ratio method, where different fragments 
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Figure 4. (Color Online) Natural logarithm of the isotopic 
yield ratio as a function of the fragment N (top) and natural 
logarithm of the isotopic yield ratio per fragment mass as 
a function of the fragment isospin asymmetry m/ (bottom). 
Left: system-to-system isoscaling and m-scaling 64 Zn+ 64 Zn, 
70 Zn+ 70 Zn. Right: isoscaling and m-scaling using different 
rris bins of the reconstructed quasi-projectile (see text). 



Figure 5. (Color Online) 




-, a and obtained from 
isobaric yield ratio (top), isoscaling (middle) and m-scaling 
(bottom) methods, applied to two different reactions ("Sys- 
ToSys") and two m s bins (Am s ). The source reconstruction 
constraints are applied to study their impact on the measured 
values. The error bars, when not visible, are comparable to 
the size of the symbols. 



considered in the calculation (see Eq. [T2")) may be affected 
differently. To investigate the effect of a good resolution 
up to Z = 18, as in Ref. [l8|], we restricted the analy- 
sis to the detectors with the highest available resolution, 
the so-called Super Telescopes. Such a restriction does 
not modify the values of W 1 , therefore it is not plotted 
in the figure. Moreover, the selection of quasi-projectile 



fragmentation events does not modify the 



trend, as 



can be seen from the figure, since stars (no data selec- 
tion), circles (Vcut), triangles (Vcut and SumZ) and full 
squares ( Vcut, SumZ and Qcut) overlap. 

We turn now to analyse the impact of m s (asterisks) 
and excitation energy (squares) constraints. The param- 
eters a and AH/T are plotted for AmJ = 0.185, while 
C S y m /T has been obtained for (m s ) = 0.183. The excita- 
tion energy was restricted to be between 3.5 and 5AMeV. 
The behavior is flat, for both a and AH/T, and presents 
less fluctuations with respect to the system-to-system val- 
ues. The average value is 0.96 ± 0.01, to be compared to 
the system-to-system average value, 0.34 ± 0.01. The 



Figure 6. Isoscaling parameters —/3 vs. a for all the Ami 
and excitation energy bins. The full line is the best fit to the 
data, while the dashed line represents a = —/3. 



difference in the values observed is due to the different 
values of A in the two cases. Indeed, the values of A, 
determined as in Eq. [5] for isoscaling between two dif- 
ferent systems and two different m s bins are 0.027 and 
0.084 respectively, which give consistent values of a/A 
(a/A = 12.6 ± 0.4 and 11.43 ± 0.12, respectively). The 
constancy of a/ A was previously shown in [40] . 
A possible excitation energy dependence of a and AH/T 
is suggested by systematically higher (even if in agree- 
ment) a and AH/T values obtained for a selected source 
excitation energy (squares), compared to the ones ob- 
tained with no energy selection (asterisks). Indeed the 
average excitation energy of the source in the latter case 
is higher (5.6±1.6AMeV) than the average excitation en- 
ergy in the selected window (4.2 ± 0.4AMeV). This trend 
is in agreement with the excitation energy dependence of 
the isoscaling parameter a/ A observed for Kr+Ni sys- 
tems in [40j. 

The selections in m s and excitation energy do not 
modify the increasing trend of a ji m extracted by the 
isobaric yield ratio method. Though the values agree 
within statistical uncertainty, the m s selection (aster- 
isks and squares) results in a systematic lower 'ji™ for 
A ^ 19. A possible excitation energy dependence may 
be observed also in jt" 1 extracted by the isobaric yield 

ratio method. Indeed for A ^ 19 3 I m values obtained 
from data selected in excitation energy (squares) are very 
slightly higher than those obtained from data not selected 
in excitation energy (asterisks). Moreover when perform- 
ing an excitation energy selection we notice that differ- 
ent mass regions are populated depending on the source 
excitation energy: higher mass fragments (A > 15) are 
produced by low excited sources, while lower mass frag- 
ments (9 < A < 15) are mainly produced by highly ex- 
cited sources. This is consistent with fragments being 
produced by different reaction mechanisms. Selection of 
heavy fragment multiplicity may weaken the dependence 
of j,™ on A. However the present statistics do not al- 
low this selection. These observations suggest that, in 



0.002 0.004 0.006 0.008 0.01 0.012 

1/A 2 

Figure 7. Quantity -1/A 2 [lnR(3, 1) - lnR(l, -1)] vs 1/A 2 
calculated from isobaric yield ratios for an average source 
isospin asymmetry ml = 0. The dashed and full lines are 
the best fits to the data with a first degree function, includ- 
ing and not including the Coulomb contribution, respectively. 

the absence of a source reconstruction, attention should 
be paid to mixing reaction mechanisms when drawing 
general conclusions. 



VI. COMPARISON OF THE THREE METHODS 

The following analysis was performed applying the de- 
scribed selections on the reconstructed quasi-projectile 
(see Sec. |V]for details). The quasi-projectiles were di- 
vided into bins based on their composition m s and their 
excitation energy per nucleon. 

The isoscaling parameters a and ft were extracted si- 
multaneously by a global fit to the yield ratios of all the 
available isotopes and, separately, by individual fits to 
the yield ratios for isotopes with a given Z. The ob- 
tained values were in good agreement. The parameter ft 
shows the same trend as a function of iV as a does as a 
function of Z, but has the opposite sign. The parame- 
ters, evaluated for all possible m s and excitation energy 
combinations, are plotted in Fig. [6l The best fit and 
the line representing a = —ft are also plotted in the fig- 
ure. The relation a = —ft appears to be approximately 
satisfied in our data [—ft — 1.11 a), thus supporting the 
equivalence of the isoscaling and the m-scaling. There- 
fore we can use Eq. [7] to calculate "l m from m-scaling, 
using a = AH/T. Small differences between a and —ft 
might be attributed to residual Coulomb effects. Indeed, 
although the ratio of the yields of a single isotope is con- 
sidered, the two sources, whose m s is fixed, might have 
different charges. We remind the reader that the SumZ 
selection constrains the source charge to be Z = 25 — 30. 

The residual Coulomb effect due to differences in 
the considered fragments charge is relevant only in the 
calculation of "Jf™ by the isobaric yield ratio method 
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and it influences mainly large mass fragments. From the 
mass formula the Coulomb energy for large Z can be 
written as 1331: 



A 



Z 2 

0.77— A 2/3 
A 2 



0.77 



(l-m) 2 A 2 / 3 (16) 



Adding this term to the free energy \t/T in Eq. [5J we see 
that a quadratic and a linear term in m are introduced 
that modify the symmetry energy coefficient and the 
external field. Also a term not dependent on m is 
introduced. The effect of introducing such corrections is 
extensively discussed in Ref. [33| . Here we concentrate 
on evaluating the Coulomb contribution, as well as the 
importance of o(m 4 ) and o(m e ) terms in Eq. [5] First of 
all we should notice that, assuming a spherical expansion, 
at low densities the Coulomb energy decreases as p 1 / 3 . 
A fit of the quantity -l/A 2 [lnR(3, 1) - lnR(l, -1)] 
allows us to estimate the fitting parameters of Eq. [S] 
and the Coulomb term. The fit results are shown in 
Fig. El The dashed and full lines are the best fit to the 
data with a first degree polynomial of l/A 2 including 
and not including the Coulomb correction, respectively. 
We remind the reader that 1 /A 2 oc m 2 . The parameters 
a' and d are fitting parameters. We found that o(m 4 ) 
are negligible, since we are dealing with relatively large 
fragments, which imply small m (but possibly large 
Coulomb). The Coulomb parameter d is about two 
orders of magnitude smaller than the l/A 2 coefficient 
a', but the term has to be included in the fitting func- 
tion to obtain a good fit (\ 2 /ndf ~ 0.8). We therefore 
took into account this correction and modified Eq. 1 131 as: 



C sym {T) 



A 

' 8 



lnR(3, 1, A)- 



-lnR(l,-l,A) -5(3,1, A)} -dA 2 l z 



(17) 



The a parameters are related to 3 j! m through the 
value of A, as given in Eq. [7] Therefore, the extracted 
value of the symmetry term depends on the choice of 
A. The isobaric yield ratio method instead provides a 
value independent of the choice of A. Three dif- 
ferent definitions of A have been examined. The Z/A 
of the emitting source has been adopted for the defini- 
tion of A source (Eq. [HJ. The average proton content of 
each fragment with mass greater than 4 has been used to 
calculate Au qu id (Eq. [5]) . Finally the average of the frag- 
ment (rnj) over the events has been used to define A{ mf ) 
(Eq. [Toll . The three different A values as a function of the 
fragment charge for a given AmJ and excitation energy 
combination are plotted in Fig. [51 By definition, A source 
(red, dashed) and A^ m/ ^ (black, full) are independent of 



Z. We observe that A, 



0.3A, 



3 , which is close 



to the value reported in Ref. [37| . The large discrepancy 
between the value of A source and the values of An qu id 
and A( TO/ ) is due to the fact that, while both Au qu id and 
A( TO ,) are related to the fragment composition, A source is 
related to the source composition. The quantity Au qu id 
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Figure 8. (Color Online) Values of A computed according to 
the different definitions of Eqs. I51 l9land !10l vs. the fragment 
charge Z for Ami = 0.185 and 3.5 < E* t /A < 5 AMeV. 



(asterisks) shows an odd-even behavior for Z < 8, which 
can be attributed to structure effects (49] . Such odd-even 
effect is in agreement with the trend reported in (49T-l5l| . 
The overall behavior shows a rather significant decrease 



as a function of Z. which will cause 



to increase as 



a function of Z. The very high value of Au qu id for Be 
is due to the fact that 8 Be, produced in the multifrag- 
mentation process, decays before reaching the detector, 
thus (A) deviates from the "true" centroid of the isotope 
distribution. An attempt was made to estimate the 8 Be 
yield from the N = Z nuclei yields, whose corrected value 
is plotted in the figure as a circle. A plot of the yield vs. 
mass number for rrif = fragments displays a power law 
behavior with Y oc A~ T [33]. The yield of 8 Be has been 
estimated from a fit of CA~ T to the data for even-even 
N = Z nuclei, for which the pairing contribution is the 
same. The value of A changes from 0.053 to 0.033. 

Figure [5] shows a comparison between the three meth- 
ods, with different choices of A, for a given excitation 
energy (3.5 < Ef/A < 5.0 AMeV - top panels) and for 
a given source asymmetry difference (AtoJ = 0.086 - a 
and d-f panels). The "* m values obtained by the iso- 
baric yield ratio method (triangles) shows an approxi- 
mately constant behavior for A < 20, while an increase 
is observed for A — 21 and 23. This is observed for 
all the Am s combinations and excitation energy win- 
dows. The observed increase of 3 — 4 units in 



c. 



for 

A from 11 to 23 is smaller than the increase observed 
in Ref. [l8[. We would like to emphasize that one of 
the main limitations of the isobaric yield ratio method 
is that it requires a nearly "ideal" isotopic identification, 
since any contamination by other isotopes to each (Z,A) 
fragment affects different isobars by a different amount. 
This effect is particularly relevant for the yield estimate 
of heavy particles. Indeed the mass resolution decreases 
with increasing fragment size (see Fig. [TJ) lowering the 
particle identification efficiency. While this effect largely 
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cancels out when the yield ratio of the same fragment 
produced by two different sources is taken (i.e. isoscal- 



ing and m-scaling), the 



estimation for the isobaric 



yield ratio method remains affected. For this reason data 
for A > 25, when available, are not shown for the iso- 
baric yield ratio method, while they are presented for 
the isoscaling and m-scaling methods. Indeed, we do not 
observe a change of behavior of a as a function of A ex- 
tracted from the isoscaling and m-scaling methods even 
in the high mass region (see Fig. [3]). 



The trends of 



r obtained from the isoscaling and 
m-scaling depend on the definition of A, as discussed 
in Eq. [3 First, we notice that Cs ^, m obtained from the 
isoscaling with the three definitions of A are not in agree- 
ment. This is due to the differences in the values of A, 
which is clearly shown in Fig. |SJ The Cs ^ m values de- 
termined using A SOMrce (crosses) and A/ m ,\ (circles) do 
not show any A dependence, while a clearly increasing 
trend is observed for A > 15 for values determined using 
^liquid (squares). This is caused by the dependence of 
^■liquid on A, as shown in Fig. [5] For A < 15 the struc- 

values. 



ture effects affecting Au qu id are reflected in 

c 

The 3 jT" values determined using A{ mf ) from the 
isoscaling (full circles) and m-scaling (open circles) show 
a generally good agreement for all the AtoJ combinations. 
We can see that the values obtained from the m-scaling 
are systematically higher than those obtained from the 
isoscaling, which could be due to the residual Coulomb 
effects noted above. The difference decreases as the exci- 
tation energy increases, as can be seen comparing panels 
a, d, e and /. 

Comparing the top panels, we observe that the val- 
ues obtained with the isoscaling and m-scaling, using 
A( m ), do not depend on the choice of Am s . Indeed, 
the weighted average values (10.19 ± 0.04, 10.10 ± 0.20, 
10.12 ± 0.16 and 9.22 ± 0.05, 8.90 ± 0.20, 9.07 ± 0.15 
for AmJ = 0.086, 0.180 and 0.279 and extracted by m- 
scaling and isoscaling, respectively) are consistent within 
2<T. The increase in the uncertainties reflects the decrease 
of available statistics as the considered m s bins are more 
apart, i.e. when we consider the tails of the m s distribu- 
tion (see Fig. |2J|. 



An excitation energy dependence of 
trapolated from panels a, d-f: Cs y rn 



c,. 



can be ex- 
decreases from 
11.1(±0.08) to 7.7(±0.04) as the source excitation en- 

de- 

cv 



ergy increases. This trend is observed for the A( m/ ) 



rived quantities, as well as for 



derived with the 



other A definitions. As already discussed, this trend 
is in agreement with the excitation energy de pen dence 
of the isoscaling parameter a/ A observed in [4fj|. An 
excitation energy dependence of ji m is also observed 
for the m-scaling extracted quantities, which vary from 
12.47(±0.08) to 8.20(±0.04) increasing the excitation en- 
ergy. In contrast, a weaker excitation energy depen- 
dence is found for the values extracted with the isobaric 
yield ratio method, which decreases from 10.0(±0.4) to 



8.8(±0.2) increasing the excitation energy. 

We now compare the values obtained with the differ- 
ent A definitions to the values extracted from the isobaric 
yield ratio, since they are independent on the choice of 
A. The 'JZ m values determined using A source (crosses) 
are lower than those determined with the isobaric yield 
ratio method (triangles) by a factor 2 to 3 independent 
of Am s and independent of the excitation energy. The 

°jt m values determined using A[ iqui d (squares) show a 
rather good agreement for A « 18 — 22 in the central 



mass region, where the Au qmd w z-v< m/; 
it does not reproduce the data for the low mass region, 
where structure effects influence the value of Au qu id, as 
pointed out in Fig. [8j In this region, in contrast to the 
high mass region, the isotopic resolution of each fragment 
is very good and presents a small uncertainty; thus, we 



A, 



Nevertheless, 



are confident on the extracted 



values. The 



values determined using A( m A, both from isoscaling and 
m-scaling, show good agreement with the isobaric yield 
ratio values for masses up to A « 20, for all the AmJ and 
excitation energy combinations. This is consistent with 
what was recently observed by R. Tripathi et al. in Ref. 
[37]]. In this work the authors show that the position 
of the free energy central minimum, when an external 
field is present (i.e. when m s ^ [l7j), is related to 
(to/) rather than to m s . This indicates that the connec- 
tion between quantities extracted from fragment yields 
and the symmetry energy coefficient has to be found in 
the composition of the emitted fragments, taking into 
account the event multiplicity. 



The generally flat behavior of 



as a function of A 



(circles and triangles in Fig. [5]) may be interpreted as the 
weak influence of secondary decay effects on the observ- 
ables. This in turn implies that fragments are determined 
in the instability region in agreement with the theoretical 
work of Dorso and Randrup [52]. This would also be the 
case in a first-order phase transition for infinite nuclei, 
where the size of a cluster is determined by its internal 
pressure which takes into account the surface tension, 
and the external pressure due to the gas [34]. For finite 
systems the formed fragments might reach the ground 
state by emitting low-energy gamma rays and possibly a 



neutron. Also the constancy of ■ 



might indicate that 



the assumption of constant volume is reasonable and we 
can identify this quantity as the symmetry energy rather 
than the enthalpy. 



VII. CONCLUSIONS 

Three methods to extract the symmetry energy coeffi- 
cient from fragment yields were compared. The isobaric 
yield ratio method removes the dependence on the frag- 
menting source by computing the difference of the yield 
ratios of properly chosen isobars produced by the same 
source, but the dependence on individual isobar detection 
efficiencies remains. The m-scaling and the isoscaling, 
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however, retain a dependence on the source character- 
istics in the difference of the external fields of the two 
sources (AH/T), while significantly reducing the depen- 
dence on isotopic detection efficiencies. The determina- 
tion of y m in this case depends on the choice of A. 
The symmetry energy coefficient to temperature ratio, 
s ri m , was experimentally evaluated function of the 
fragment mass with three different methods. The effect 
of the source reconstruction on the s ^ m values was anal- 
ysed, showing that Air angular coverage is useful in ex- 
tracting information from the isobaric yield ratios, since a 
limited coverage impacts Cs j! m values. Improved isoscal- 
ing and m-scaling were observed when selecting quasi- 
projectile m s bins for n-rich and n-poor systems, while 

isobaric-yield-ratio-extracted values were only slightly af- 

c 

fected. A decrease in ™ m was observed with increasing 
excitation energy. 

The isoscaling parameters a and j3 were extracted as 
a function of Z for all m s and excitation energy combi- 
nations. Our data show the equivalence of the isoscaling 
and the m-scaling, since the relation a — —(3 is approxi- 



mately satisfied. 

The y"* values extracted from the a isoscaling pa- 
rameter using A( m ) are in good agreement with the iso- 

baric yield ratio ^ m values. It indicates that the con- 

c 

nection between a and ^ m has to be found in the av- 
erage fragment isotopic asymmetry. The extracted sym- 
metry term shows a generally flat trend as a function of 
A for the mass region A = 10 — 20 with the best mass 
resolution, independent of the method used. This may 
be consistent with the lack of secondary de-excitation ef- 
fects in our data, within our isotopic resolution, as well 
as with the assumption of a freeze-out volume where the 
disassembly occurs. 
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Figure 9. (Color Online) Values of s j!"' obtained with the three different methods. Each panel corresponds to the noted Am s 

and the noted window in E* /A. Three different choices of A have been used to compute "j! m from the isoscaling parameter 
a. 



